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Abstract 

We study localized traveling waves and chaotic states in strongly nonlinear one-dimensional 
Hamiltonian lattices. We show that the solitary waves are super-exponentially localized, and present 
an accurate numerical method allowing to find them for an arbitrary nonlinearity index. Com- 
pactons evolve from rather general initially localized perturbations and collide nearly elastically, 
nevertheless on a long time scale for finite lattices an extensive chaotic state is generally observed. 
Because of the system's scaling, these dynamical properties are valid for any energy. 
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I. INTRODUCTION 



Hamiltonian lattices are one of the simplest objects in nonlinear physics, nevertheless they 
still elude full understanding. Already the first attempt to understand nonlinear effects ended 
with the Fermi-Pasta-Ulam puzzle, which is still not fully resolved (see e.g. the focus issue 
on "The "Fermi-Pasta-Ulam" problem - the first 50 years" in Jj); another remarkable feature 
found only recently is the existence of localized breathers [2||. Quite often nonlinear effects 
in lattices can be treated perturbatively, leading to well-established concepts of phonon 
interaction and weak turbulence. Beyond a perturbative account of a weak nonlinearity 
one is encountered with genuine nonlinear phenomena, like solitons and chaos. The level of 
nonlinearity usually grows with the energy, allowing one to follow a transition from linear 
to nonlinear regimes by pumping more energy in the lattice. 

In this paper we study strongly nonlinear Hamiltonian lattices that do not possess linear 
terms. We restrict our attention to the simplest one-dimensional case where particles inter- 
act nonlinearly and no on-site potential is present. We choose the interaction potential in the 
simplest power form, thus the lattice is characterized by a single parameter - the nonlinear- 
ity index. The equations of motion obey scaling, what means that the dynamical properties 
are the same for all energies - only the time scale changes. Lattices of this type attracted 
large attention recently, in particular due to a prominent example - the Hertz lattice that 
describes elastically interacting hard balls, it has the nonlinearity index 3/2 . We 

focus our study on the interplay of solitary waves and chaos in such lattices. Some 25 years 

n p n Q 

ago Nesterenko [3|, lZ|, 18|, 19|| has described a compact traveling wave solution in the Hertz 
lattice, which can be understood as a compacton. Compactons have been introduced, in a 



mathematically rigor form, by Rosenau and coworkers [10|, LUfl for a class of nonlinear PDEs 
with nonlinear dispersion. Compactons can be analytically found if one approximates the 
lattice equations with nonlinear PDEs, but less is known on the genuine lattice solutions. 
Below, in Section HH BI we present a numerical procedure for determining traveling waves for 
an arbitrary nonlinearity index, and compare these solutions with those of the approximated 
PDEs (Section IHI All . Furthermore, we show that compactons naturally appear from local- 
ized initial perturbations and relatively robustly survive collisions, but nevertheless evolve 
to chaos on a long time scale in finite lattices (Section [TV]). Properties of chaos are studied 
in Section [Vl We demonstrate extensivity of the chaotic state by calculating the Lyapunov 
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spectrum, and study the dependence of the Lyapunov exponents on the nonlinearity index. 
Some open questions are discussed in the concluding Section |VTl 



II. THE MODEL 



Our basic model is a family of lattice Hamiltonian systems 
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which are parameterized by one real parameter - the nonlinearity index n. Below we assume 
that n > 1. The case n = \ corresponds to a linear lattice. Another interesting case is 
n = 3/2. Such a nonlinearity appears, according to the Hertz law, at the compression 
in a chain of elastic hard balls. For the realistic system of balls, however, the potential 
has the form like in ^ only for g^+i — qk < 0, for qk+i — > no attracting force is 
acting. A simplified realization of such a system is the toy "Newton's cradle", which possesses 
the same Hertzian interaction law. However the standard "Newton's cradle" consists of 
a few balls (typically 5), which are not enough for the formation of stationary traveling 
waves. Furthermore, slight intervals between adiacent beads are not excluded, contrary to 

HQ 

experiments [4, iZfl where great care is taken to let the beads in effective contact. For different 
ispects of the Hertz chain see refs. 0, ^, [l2, 13, 1^, 15, 16, 17, 1^, 19, 20, 21|, a review article 
22|, and references therein. Contrary to this, in our model ([T]) we assume both repulsive 

and attracting forces. 

Note that the potential in ([T]) is generally non-smooth, except for cases n = 1,3,5, .. .. 
Although the dynamics can be easily studied in non-smooth situations as well, we will mainly 
focus below on the simplest smooth non-trivial case n = 3. 

The lattice equation of motion read 



Qk 



\qk+i - gfcrsign(gfc+i - Qk) 



\Qk 



gfc_i|''sign(gfc - g^. 



(2) 



Since on the right hand side of (l2|) only differences enter, it is convenient to introduce the 
difference coordinates Qk = qk+i — Qk- Then the equations of motion are transformed to 



Qk = IQk+if sign{Qk+i) - 2 IQkl"' sign{Qk) + IQfc-iT sign(Q 



k-li 



(3) 



Note, that a solitary wave in the variables Qk corresponds to a kink (shock-like wave) in the 
variables Qk- 



Conservation laws. The equations of motion possess two conservation laws: the energy 
and the total momentum, the latter can be trivially set to zero by transforming into a moving 
reference frame. 

Scaling. As mentioned in Refs. , the lattice ^ has remarkable scaling properties, 

due to homogeneity of the interaction energy. It is easy to check, that the Hamiltonian can 
be rescaled according to 

q = aq, p = a~p, H = H, t = a~t. (4) 

Note, that this scaling involves only the amplitude and the characteristic time of the so- 
lutions: by decreasing the amplitude one obtains new solutions having the same spatial 
structure but evolving slower. We will see that this property has direct consequences for the 
properties of traveling waves and of chaos. 



III. TRAVELING SOLITARY WAVES 



In this section localized traveling waves are investigated, first in a quasicontinuous approx- 
imation (QCA) and then via numerical solution of the lattice equations. A mathematically 
rigor proof of the existence of solitary waves in Hamiltonian lattices of type ^ has been 
given in Refs. 



A. Quasicontinuous approximation 

Here, we represent the solution of the lattice equations ((31) as a function of two continuous 
variables Q{x,t). We are seeking for solitary waves which do not change their sign, for 
definiteness we consider Q > (this consideration is therefore suitable for lattices where 
the nonlinearity index is different for positive and negative displacements Q, e.g. for the 
Hertz lattice of elastic balls) . We present two approaches to find a continuous version of the 
lattice. In the first one, we approximate the differences between two displacements Q, while 
in the second one the displacement q at each lattice site is expanded directly. 
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1. Expansion of differences 



Here we look for a direct quasicontinuous approximation of Eq. (I3|). Expanding the 
difference coordinates Qk up to the fourth order we obtain 



where h is the spatial difference between two lattice sites and the subscripts denote differen- 
tiation with respect to x. Inserting ((51) into ((31) and setting h = 1 one arrives at the partial 
differential equation 

[Q{x. t)]u = mx, t)i^ + ^ mx, t)]_ . (6) 

Equation ((61) belongs to a class of strongly nonlinear PDEs, because the dispersion term with 
the fourth derivative is nonlinear. The equation does not possess linear wave solutions (this 
situation has been called "sonic vacuum" by V. Nesterenko), but it has nontrivial nonlinear 
ones. In this way it is very similar to a family of strongly nonlinear generalizations of the 
Korteveg-de Vries equation, studied in {l^ and can be considered as a strongly nonlinear 
version of the Boussinesq equation 

Now we seek for traveling wave solutions of ([6]), by virtue of the ansatz 

Qix,t) = Q{x- \t) = Q{s). (7) 

Then ^ reduces to the ODE 

>^'Qss=[Q\s + ^mssss- (8) 

Furthermore, we assume that the solution tends to zero as s ^ ±oo, thus after integrating 
twice we obtain 

A'Q = ^" + 1^ [Q^'ls ■ (9) 
This equation also appears in the traveling waves ansatz for the (n, n)-equation in (lo| . 
Eq. ^ can be solved for an arbitrary power n by 

g(s) = |ArAicos™(5is) (10a) 

with ^ 

m= -, Ai=i^—] , B, = V3 . 10b 

n — 1 \ 2n I n 
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2. Expansion of displacements 



Another type of quasi-continuum can be obtained if we approximate Eq. Now the 
displacement q at each lattice site is written as a continuous variable, what for the same 
order of the spatial derivative as in ([5]) gives: 

qk±i = q±hq^ + —q^x ± -prqxxx + t^I^Ixxxx (H) 

Inserting this expansion into the equations of motion ((21) and collecting all terms up to order 
of K^^^ yields 

|9]« = ft"*' Kl. + ^ - . (12) 

TMs equation is the long wave app.ci^a.ion of Neste.ento fl 0. To co.pa. « with 

Eq. ([6]), we differentiate (fT2ll with respect to x, define Q = hq^ and set = 1: 

= [Q%x + ^ - '^^^^^^[Q'^-'QIU) . (13) 

One can see that there is an additional term in (fT3ll compared to ((61). This is not so 
much surprising, as these two quasi-continuous approximations correspond to expansions at 
the different positions of the original lattice, this effect is well-known for approximations 
of Hamiltonian lattices with PDEs 3|- Because in the problem we do not have a small 
parameter (the lattice spacing h = 1 \s not small compared to the wave length), none of 
the equations (l6ll. (fT3l) can be expected to be exact in some asymptotic sense. Instead, one 
has to justify them by comparing the solutions with those of the full lattice problem, see 
Section IIII Bl below. 

To find traveling waves in the direct expansion we use again the ansatz ((71) Q{x,t) = 
Q{x — Xt) = Q{s). Inserting this ansatz and integrating twice yields then an analogon to ([9]) 

= 0" + ^mss - (14) 

One partial solution of this ordinary differential equation can also be written as 

Q{s) = \XrA2 cos"" {B2S) (15a) 



but with different constants A2 and B2 (cf. |l9l. l26l|) 



m = -, A2=[-—) , B2 = \Q \ ' 15b 

n — 1 \l + n J y n[n + 1) 



The solutions (iTOl) and (1151) do not satisfy boundary conditions, moreover, they intersect 
with another, trivial solution of (l8|) Q = 0. Remarkably, because of the degeneracy of Eq. ((HI) 
and Eq. f fT4l ) at zero, one can merge periodic solutions (fTOl) and f fTS! ) with the trivial solution 



Q = (see a detailed discussion in 



Q{s) 
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III): 



|A|'^A cos"" {Bis) 



s < 



else. 



2Bi 



(16) 



with i = 1,2. This gives a compacton - a solitary wave with a compact support - according 



to definition 







lll |. For other, non-solitary solutions of ( fT3l ). see e.g. [9|, l26||. Note that 
due to the symmetries x — x, Q —Q, solitary waves with both signs of velocity A and 
of amplitude A are the solutions. 

It is important to check the validity of solution (flGll by substituting it back to ((HI) or (fT3l) . 
Then no terms are singular for the case m > 2 only, i.e. for n < 2. Thus, the constructed 
compacton solution (fTBl) is, strictly speaking, not valid for strong nonlinearities n> 2. This 
conclusion is, however, only of small relevance for the original lattice problem. Indeed, the 
PDE ([6]) or f |T2l ) is only an approximation of the lattice problem: because the spatial extent 
of solution dHj) is finite, there is no small parameter allowing us to break expansion (l5l) or 
(fTTj) somewhere. Just breaking it after the fourth derivative is arbitrary and can be justified 
only by the fact, that in this approximation one indeed finds reasonable solutions at least for 
some values of n. A real justification can come only from a comparison with the solutions 
of the lattice equations, to be discussed in the next subsection. And there we will see that 
the solution can be found both for weak and strong nonlinearities n > 2. 



B. Traveling waves in the lattice 

In the lattice, the traveling wave ansatz reads Qk(t) = Q{k — \t) = Q{s). Inserting this 
ansatz into the lattice equations ([3]) yields 

\^Q"{s) = Q^{s - 1) - 2Q"(s) + Q"(s + 1). (17) 



We employ now the scaling (Ilj) and set A = 1. As demonstrated in 
differential equation can be equivalently written as an integral equation 



27|, this advanced-delay 



Q{s)= f (l-k-^|)g"(Ode- (18) 

Js-l 



1.6 




(e) (f) 

FIG. 1: The traveling waves obtained from (fTOl) for various powers n. Markers show the wave on 
the lattice, dotted lines show the corresponding solutions of the quasi-continuous approximation 
(ID) and dashed lines show the solutions of the QCA (fT3ll . Left column: normal scale, right column: 
logarithmic scale. (a,b) n = 3/2; (c,d) n = 3, (e,f) n = 11. Note that the width w of the compacton 
decreases as n increases. 
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[One can easily check the equivalence by differentiating (fTSll twice.] We can now, following 



the approach of V. Petviashvili 



28 



29l |. construct an iterative numerical scheme to solve the 



integral equation ( fTSl) . Starting with some initial guess Qq, one constructs the next iteration 
via 



Qi+l 



\\Q. 

+1 



Q* and 



(19a) 
(19b) 



-1 



(practically, we used the Li-norm for || ■ ||). We have used a = that ensured convergence 
of the iterative scheme. The integral in (fT9l) was numerically approximated by virtue of a 
4th-order Lagrangian integration scheme [30|. In Fig. [T] the traveling waves for various 



powers n are shown. Using the logarithmic scale one clearly recognizes the compact nature 
of the waves. 

In Fig. [2] we show the dependencies of the total energy E, the solution Li-norm Nl^, 
and the amplitude Qmax of the found waves on the nonlinearity index n, for a fixed wave 
velocity of A = 1. Remarkably, the effective width Ni^/Qmax decreases with increasing 
nonlinearity index and it seems that the profile of the compacton converges to a triangular 



20||, where the dependence of the 



shape as n — i> oo. A similar result has been obtained in 
pulse velocity on the nonlinearity index has been analyzed for large n in the binary collision 
approximation. 



C. Estimation of the tails 



It is clear from the integral form (ITSll . that the solution cannot have a compact support. 
In this section we estimate the decay of the tails. We start with (ITSll and substitute Q{s) = 

s+l 

Qis) = I (l-|s-e|)e-"^«)de (20) 

We consider the tail for large s > 0, if we assume a rapid decay of Q(s), then the integrand 
in ([20l) has a sharp maximum at s — 1. Thus we can approximate the integral using the 
Laplace method. At the maximum we expand /(^) into a Taylor series around s — 1, keeping 
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n 



FIG. 2: The dependency of the amplitude Qmax, the energy E and the Li-norm Nl^^ of a com- 
pacton on the nonhnearity index n. In this plots A = 1. For comparison, the curves from the 
quasicontinuous approximation are shown with dotted lines for the Eq. jH]) and with dashed lines 
for Eq. 

only the leading first-order term: 

s+l 

Q{s) ^ / (1 - 1^ - i\) exp { - nf{s - 1) - nf\s - - (s - l))]di (21) 



We shift the integration range 



Q{s) ^ e 



(22) 



where we also replace the decreasing part of the kernel with ^. Since this integrand decreases 
very fast, we can set the upper bound of the integration to infinity, then by partial integration 
we obtain 



Q{s) = e 



-nf{s-l) 



[nf'is - 1 



Taking the logarithm of this equation yields 



- f{s) = -nf{s - 1) - 21og[n/'(s - 1)] 



(23) 



(24) 
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FIG. 3: (a) The tails of the compactons in double logarithmic scale, (b) Comparison of the estimate 
([271) with the compactons obtained from (fTOl) . 



Since we expect that f{s) is a rapidly growing function of s, we can neglect the logarithmic 
term and obtain 

fis) = nf{s - 1) (25) 

This equation is solved by 

/(s) = Cn' = Ce^°g(")^ (26) 
where C is an arbitrary constant. Finally we obtain that the tail decays super-exponentially: 



Q{s) = e~^(^) ^ e-^"' = exp 



C exp ( log(n)s) 



(27) 



This expression was first obtained by Chatterjee 
delayed equation (flTll . 



5| using a direct expansion of the advanced- 



In Fig. 3(a) we show the tails of the compactons for various values of n and in 



Fig. 3(b) we compare the estimated decay rate (l27l ) with compactons obtained numeri- 



cally from the traveling wave scheme f |T9l l. To obtain the double logarithmic decay rate 
P = d\og{\ \og{Q{s))\)/ds, we first compute log(| log(Q(s))|) and then the derivate is calcu- 



lated using a spline smoothing scheme 



3l| . To suppress small oscillations of the tails we 



average the numerical obtained derivative in the last 1/6 of the compacton domain. The 



numerical value of P is shown in Fig. 3(b) Both coincide very well. 
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IV. EVOLUTION AND COLLISIONS OF COMPACTONS 



A. Appearance of compactons from localized initial conditions 

The compact solitary waves constructed in the previous section are of relevance only if 
they evolve from rather general, physically realizable initial conditions. For an experimental 
significance (see 0, Q| for experiments with Hertz beads) it is furthermore important, that 
the emerging compact waves establish on relatively short distances, otherwise dissipation 
(which has not been considered here) will suppress their formation. We illustrate this in 
Figs. [4f5l There we report on a numerical solution of the lattice equations ((21) on a finite 
lattice of length N = 128 (so that at the boundaries g'l = \q2 — gi|"^sign(g2 — ^i) and 

= —\qN — qN-i\"'siga{qN — qN-i) holds). One of the quantities we report is the local 
energy at site k defined as 

^fc = y + ^^^^ - Qk^' + \qk - qk-ir') • (28) 

As an initial condition we have chosen a kink in the variables g^: = (n+l)^/*^"+^) for /c > 64 
and qk = elsewhere. This profile has unit energy, it corresponds to the localized initial 
condition in the variable Q: Qk = Sk,m ■ in' + l)i/("+i). The evolution of different variables is 
shown in Fig. [H From the initial pulse of Q, a series of compactons with alternating signs is 
emitted in both directions. The amplitude of the perturbation near the initially seeded site 
decreases and correspondingly increases a characteristic time of the evolution. We expect 
that at large times, compactons with small amplitudes will continue to detach. In Fig. [5] we 
show the evolution from the initial step for different nonlinearities n = 1.5, 3, 10. The plots 
look very similar and compactons are emitted in every case. The number of the emitted 
compactons and their amplitudes depend on the nonlinearity index. 

In our next numerical experiment we studied the emergence of compactons not from a 



sharp step in the coordinates g^, but from localized random initial conditions. In Fig. 6(a) 
we show a typical evolution in a lattice of length = 512 (with nonlinearity index n = 3) 
resulting from random initial conditions qk in the small region N/2 — 5 < k < N/2 + 5 around 
the center of the lattice. In this region the coordinates qk have been chosen as independent 
random numbers, identically and symmetrically uniformly distributed around zero, while 
Pk{0) = 0. Furthermore, the energy of the lattice was set to = 1 by rescaling. In a 



particular realization of Fig. 6(a), at the initial state two compactons emerge to the right 
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FIG. 4: Evolution from an initial step for the nonlinearity index n = 3, the lattice length is = 128 
and open boundary conditions (hence q'l = {q2 — qi)"^ and q'N = —{qN — QN-i)'^) are used. The 
initial conditions are qk{t = 0) = (n + for k > 64 and else, initial momenta are zero. 

Different plots show different quantities of the lattice: (a) the coordinates q^; (b) the energy 
defined in ([28]l : (c) the difference coordinates Qk = qk+i — Qk at time t = 80, the initial state at 
t = is shown here as the dashed line and (d) the difference momenta Pk = Pk+i — Pfc at t = 80. 
The compactons originating from this initial state are clearly separated near the borders of the 
chain, those in the middle part are still overlap-ed. 

and four compactons to the left. In the center of the lattice a chaotic region establishes 
and slowly spreads over the lattice, possibly emitting more compactons on a longer time 



scale. In Fig. 6(b) v^e perform a statistical analysis of this setup by showing the energy 
distribution of compactons emitted from localized random initial conditions as described 
above. This distribution was obtained from 60000 simulations, in each simulation the energy 
of the emitted compactons have been determined and counted. The functional form of the 
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FIG. 5: Evolution of an initial step (like in Fig. |4| for various nonlinearity indices (a) n = 1.5, (b) 
n = 3 and (c) n = 10. The initial state is shown as dashed line, the solid line is the state at time 
t = 80. 




(a) (b) 



FIG. 6: (a) Compactons emerging from localized random initial conditions. The nonlinearity index 
is n = 3. The gray scale corresponds to the energy ([28]l of the lattice site, (b) Energy distribution 
of the compactons emitted from localized random initial conditions. The statistics was obtained 
from 60000 simulations; in each simulation the lattice was integrated to the time T = 1000 and 
the energy distribution of the compactons emerged to the right (black circles) and the left (crosses) 
have been determined. The distributions obeys in very good approximation P{E) ~ E~<^^°siE)-b ^ 
with a = 0.57 and b = 5.47. 

distribution obeys in very good approximation P{E) ~ ^-«iog(^^)-*^ with a = 0.57 and 
b = 5.47. 
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B. Collisions of compactons 

As we have demonstrated above, compactons naturally appear from rather general initial 
conditions. To characterize their stability during the evolution, we study their stability to 
the collisions. This study is not complete but only illustrative, as in Fig. [7] we exemplify 
different cases of collision of two compactons in a lattice with n = 3. These six setups 
present all possible scenarios of two compactons: (i) two colliding compactons with the 
same amplitudes, (ii) two compactons with different amplitudes moving toward each other, 
and (iii) two compactons with different amplitudes moving in the same direction and passing 
each other. Each of these three cases has two sub-cases, because the amplitudes can have 
the same or different sign. It should be mentioned, that these six collisions do not represent 
the complete picture of all collision, moreover we have not varied parameters such as the 
distance between two colliding compactons or their amplitudes. 

In all the cases presented, the initial compactons survive the collision: they are not 
destroyed although they do not come out of the collision unchanged. In all cases the collision 
is non-elastic, some small perturbations (that presumably on a very long time scale may 
evolve into small-amplitude compactons) appear. 

Because of this non-elasticity, on a finite lattice after multiple collisions initial compactons 
get destroyed and a chaotic state appears in the lattice, as illustrated in Fig. [3 There we 
show the evolution of the two compactons with the same amplitude and sign of the amplitude 
for three different nonlinearities: n = 3, n = 9/2 and n = 11. In the first two cases the 
chaotic state establishes relatively fast. In the third simulation with n = 11 the situation is 
different. Here the chaotic state does not appear even on a very long time scale. We run 
the simulation for very long times up to T = 2 ■ 10^, but could not observe the development 
of a chaotic state. We have checked this phenomenon also for higher values of n with the 
same result. Presumably, these initial conditions lie on a stable quasiperiodic orbit or are 
extremely close to a such one. 

V. CHAOS IN A FINITE LATTICE 

As demonstrated above, in a finite lattice general initial conditions evolve into a chaotic 
state. For characterization of chaos we use Lyapunov exponents. The chaotic state of the 
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FIG. 7: Collisions of compactons in the Hamiltonian lattice with n = 3, shown are difference 
coordinates Qk- We have considered 6 different collision scenarios, in each plot the lower panel is 
the initial configuration of the lattice, the middle panel is the state of the lattice at some time during 
the maximal overlap, and the upper panel shows the lattice past the collision, (a) Compactons of 
equal energy having opposite amplitudes and velocities; (b) compactons of equal energies and 
amplitudes but opposite velocities; (c) compactons of different energies having amplitudes and 
velocities of opposite signs; (d) compactons of different energies having amplitudes of the same and 
velocities of opposite signs; (e) compactons of different energies having velocities and amplitudes of 
the same sign; (f) compactons of different energJ^ having velocities of the same sign and amplitudes 
of opposite signs. 




FIG. 8: Collisions of compactons and emergence of chaos after multiple collisions. Different plots 
show different nonlinearity indices (a) n = 3, (b) n = 9/2 and (c) n = 11. Time increases from left 
to right and the difference coordinates Qk are shown in gray scale. Remarkably the elasticity of the 
collision increases with increasing nonlinearity index n, so that practically no irregularity appears 
at n > 10. 



lattice has also been characterized in 



22 



32l | by the means of the velocity distribution of 
the lattice site. It has been found that the lattice possesses a quasi-non-equilibrium phase, 
characterized by a Boltzmann-like velocity distribution but without energy equipartition. 
First, we check that chaos in the lattice is extensive, i.e. the Lyapunov exponents form a 



spectrum when the system size becomes large (Fig. 9(a)). This property allows us to extend 
the calculations of finite lattices to the thermodynamic limit. Note, that due to the two 
conservation laws, four Lyapunov exponents vanish; we have not found any more vanishing 
exponents, indicating the absence of further hidden conserved quantities. 

For a lattice of length = 16 the dependence of the Lyapunov exponents on the non- 



linearity is shown in Fig. 9(b) For a fixed total energy (we have set H = N = 16 in these 
calculations) the Lyapunov exponents grow with the nonlinearity index. The plot presented 



in Fig. 9(c) indicates that Xmax logn, although we did not consider very high nonlinearity 
indices to make a definite conclusion on the asymptotics for large n. 

We stress here that because of the scaling of the strongly nonlinear lattices under con- 
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index k/(N-2) index k/(N-2) n 

(a) (b) (c) 

FIG. 9: Lyapunov exponents of the Hamiltonian (a) The Lyapunov spectra for one fixed 

nonhnearity index n = 3 and different values of lattice length N = 16, 32, 64, 128. The index axis 
is normalized to 1. (b) The Lyapunov spectra Xj for various values of the nonlinearity index n 
(from bottom to top: n = 1.5, 2, 2.5, 3, 3.5, 4, 6, 8) and fixed lattice length N = 16. Larger values 
of n produce stronger chaos that smaller ones, (c) The largest Lyapunov exponent Ai for different 
values of n. The horizontal axis is logarithmic, thus one can see that roughly Aq ~ const • log(n). 

sideration, chaos is observed for arbitrary small energies - only the Lyapunov exponents 
decrease accordingly. 



VI. CONCLUSION 



In this paper we have studied strongly nonlinear Hamiltonian lattices, with a focus on 
compact traveling waves and on chaos. We have presented an accurate numerical scheme 
allowing one to find solitary waves. Moreover, from the integral form representation used 
one easily derives the super-exponential form of the tails. In this way we have confirmed this 
remarkable result by Chatterjee by another analytical method and by accurate numerical 
analysis. The constructed compactons were then studied via direct numerical simulations of 
the lattice. Their collisions are nearly elastic, but the small non-elastic components on a long 
time scale destroy the localized waves and result in a chaotic state. Chaos appears to be a 
general statistically stationary state in finite lattices, with a spectrum of Lyapunov exponents 
where the largest one grows roughly proportional to the logarithm of the nonlinearity index. 

We would like to mention here also several aspects that deserve further investigations. 
Recently, a problem of heat transport in one-dimensional lattices have attracted a large 



attention 



33|, here the properties of strongly nonlinear lattices may differ from those pos- 



sessing linear waves. Also a quantization of these lattices seems to be a non-trivial task, as 
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there are no linear phonons to start with. Finally, the Anderson localization property of dis- 
ordered lattices has been recently intensively discussed for nonlinear systems. For strongly 
nonlinear lattices the problem has to be attacked separately, as here one cannot rely on the 
spectral properties of a linear disordered system. 
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